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Abstract 

Here we study a class of second-order nonautonomous differential equations, and 
the corresponding planar and spatial systems, from the point of view of fractal 
geometry. The fractal oscillatority of solutions at infinity is measured by oscillatory 
and phase dimensions. The oscillatory dimension is defined as the box dimension 
of the reflected solution near the origin, while the phase dimension is defined as 
the box dimension of a trajectory of the corresponding planar system in the phase 
plane. Using the phase dimension of the second-order equation we compute the 
box dimension of a spiral trajectory of the spatial system, lying in Lipschitzian or 
Hblderian surfaces. Also, the box dimension of a trajectory of the reduced normal 
form with one eigenvalue equals to zero, and a pair of pure imaginary eigenvalues 
has been computed at the moment of Hopf bifurcation. 



1 Introduction and motivation 

According to the idea of qualitative theory of differential equations, oscilla- 
tions of a class of second-order differential equations have been considered by 
phase plane analysis, in Pasic, Zubrinic and Zupanovic [Hj. The novelty was 
a fractal approach, connecting the box dimension of the graph of a solution 
and the box dimension of a trajectory in the phase plane. Oscillatory and 
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phase dimensions for a class of second-order differential equations have been 
introduced. The notion of fractal oscillatority near a point for real functions of 
a real variable has been introduced in Pasic [TT] , and studied for second-order 
linear differential equations of Euler type, as the basic model. The study was 
based on the fact that the Euler equation has chirp-like oscillatory solutions. 
This fractal approach to the standard Bessel equation and a generalized Bessel 
equation, can be found in Pasic and Tanaka [T2"iri3] . and also in [9]. 

On the other hand, in Zubrinic and Zupanovic [19.22J, the box dimension of 
spiral trajectories of a system with pure imaginary eigenvalues, near singular 
points and limit cycles, has been studied using normal forms, and the Poincare 
map. These results are based on the fact that these spiral trajectories are of 
power type in polar coordinates. Also, it is shown that the box dimension is 
sensitive with respect to bifurcations, e.g., it jumps from the trivial value 1 
to the value 4/3 when the Hopf bifurcation occurs. Degenerate Hopf bifur- 
cation or Hopf-Takens bifurcation can reach even larger box dimension of a 
trajectory near a singular point. This value is related to the multiplicity of 
the singular point surrounded with spiral trajectories. This phenomenon has 
been discovered for discrete systems in Elezovic, Zupanovic and Zubrinic [1] 
concerning saddle-node and period-doubling bifurcations, and generalized in 
Horvat Dmitrovic [5]. Also, in Mardesic, Resman and Zupanovic [10] there 
are results about multiplicity of the Poincare map near a weak focus, limit 
cycle, and saddle-loop, obtained using the box dimension. Isochronicity of a 
focus has been characterized by box dimension in Li and Wu [6] . Formal nor- 
mal forms for parabolic diffeomorphisms have been characterized by fractal 
invariants of the e— neighborhood of a discrete orbit in Resman [15]. All these 
results are related to the 16th Hilbert problem. 

This work is a part of our research, which has been undertaken in order to un- 
derstand relation between the graph of certain type of oscillatory function, and 
the corresponding spiral curve in the phase plane. We believed that chirp-like 
oscillations defined by A(r) = r a sinl/r "correspond" to spiral oscillations 
r = (f~ a in the phase plane, in polar coordinates. The relation between these 
two objects has been established in [8], introducing a new notion of the wavy 
spiral. Applications include two directions. Roughly speaking, we consider spi- 
rals generated by chirps, and chirps generated by spirals. If we know behavior 
in the phase plane, we can obtain the behavior of the corresponding graph, 
and vice versa. As examples we may consider weak foci of planar autonomous 
systems, including the Lienard equation, because all these singularities are 
of spiral power type r = ip~ a , a G (0,1), see [T4"fT9"] . As an application of 
the converse direction, from a chirp to the spiral, we were looking for the 
second-order equation with an oscillatory solution having chirp-like behavior. 
The Bessel equation of order v is a nice example of such behavior, see [9]. 
Whereas the Bessel equation is a second-order nonautonomous equation, we 
interpret the equation as a system in IR 3 , using t — > oo, instead of the standard 
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approach with a variable near the origin. The system studied in this article, 
see (|2j) , coincides with the Bessel system for p(t) = t~ a , a = v = 1/2. We 
classify trajectories of the system with respect to their geometrical and fractal 
properties. 

Also, the results about spiral trajectories in IR 3 , from Zubrinic and Zupanovic 
[2T)f2"Tj are extended to the systems where the Hopf bifurcation occurs. The box 
dimension of a trajectory of the reduced normal form with one zero eigenvalue, 
and a pair of pure imaginary eigenvalues, has been computed at the moment 
of the Hopf bifurcation. Essentially, the Hopf bifurcation studied here is a 
planar bifurcation, but the third equation affects the box dimension of the 
corresponding trajectory in the space. We show that in 3-dimensional space, 
the Hopf bifurcation occurs with the box dimension of a spiral trajectory larger 
than 4/3, which is the value of the standard Hopf bifurcation. 



Our intention is to understand a fractal connection between oscillatority of 
solutions of differential equations and oscillatority of their trajectories in the 
phase space. Our work is mostly motivated by two nice formulas from the 
monograph of C. Tricot [181 P- 121]. He computed the box dimension for a 
class of chirps and for a class of spirals of power type in polar coordinates. We 
are looking for a model to apply these formulas, and also to show that chirps 
and spirals are a different manifestations of the same phenomenon. Here we 
study, as a model a class of second-order nonautonomous equations, exhibiting 
both chirp and spiral behavior 
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where p is a function of class C 2 . The explicit solution is x(t) 
C 2 p(t) cost. If z — l/(t — C 3 ) we get the cubic system 
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In this article we compute the box dimension of a spiral trajectory of the 
system ^ exploiting the dimension of (a, 1)— chirp X(r) = r a sinl/r, a £ 
(0, 1), for r > small. Using change of variable for time variable r 1— > r -1 , the 
infinity is mapped to the origin, and such reflected solution of [T] with respect 
to time is called the reflected solution. We use notation t = r -1 . If the function 
p(t) in ([!]) is "similar" to the function t~ a , then the reflected solution of x(t) 
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X(t) = C l p{ 1 ) sin r- 1 + C 2 p(-) cosr- 1 
r r 

1 CI 

C\ + C| p(-) sin(T _1 + arctan -^), r G (0, — ], 
r Ci t 

is an (a, l)-chirp-like function near the origin, see [7]. We obtain results con- 
necting functions "similar" to (a, 1)— chirps, and spirals "similar" to r = ^p~ a , 
a G (0,1), in the phase plane. Applications include nonautonomous planar 
systems, so we introduce the third variable z depending on the time t. Fur- 
thermore, if we take z = t~ 7 , 7 > 0, in (JTl), the box dimension of a tra- 
jectory depends on 7. For some values of 7 trajectory in IR 3 is obtained as 
bi-Lipschitzian image of the spiral from the phase plane, which does not affect 
the box dimension. For other values, trajectory lies in the Holderian surface, 
affecting the box dimension. The Holderian surface has an infinite derivative 
at the origin, which is the point of accumulation of the spiral. Spirals of the 
Holderian type have the "tornado shape" with a small bottom and wide top. 

It is interesting to notice that our results about the box dimension of planar 
trajectories of a system with pure imaginary eigenvalues, show that the box 
dimension of any trajectory depends on the exponents of the system. In IR 3 
we have already found an example, see [20||2T] . where dimension depends on 
the coefficients of the systems, which will be the case in @. See for the 
computation of the box dimension of the system 

r = airz, ip = 1, z = b 2 z 2 , (3) 



in cylindrical coordinates. If ai/b 2 G (0, 1] then any spiral trajectory r of ^ 
has the box dimension dim^ T = fq—y^ near the origin. 



2 Definitions 



Let us introduce some definitions and notation. For A C M. N bounded we define 
£ -neighborhood of A as: A £ := {y G WL N : d(y, A) < e}. By lower s- dimensional 
Minkowski content of A, s > we mean 



M'JA) :=liminf 



and analogously for the upper s -dimensional Minkowski content A4* S (A). The 
lower and upper box dimensions of A are 

dim p A := inf{s > : M&A) = 0} 
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and analogously dim^A := inf{s > : Ai* s (A) = 0}. If these two values 
coincide, we call it simply the box dimension of A, and denote by dime A 
It will be our situation. If < Mi{A) < M* d (A) < oo for some d, then we 
say that A is Minkowski nondegenerate. In this case obviously d = dim^ A. In 
the case when lower or upper d- dimensional Minkowski contents of A are or 
oo, where d = dime A, we say that A is degenerate. For more details on these 
definitions see e.g. Falconer [2], and [T9] . 

Let x : [to, oo) — >■ R, to > 0, be a continuous function. We say that x is an 
oscillatory function near t = oo if there exists a sequence tk \ oo such that 
x(tk) = 0, and functions x\( tk) t k+1 ) intermittently change sign for k G N. 

Let u : (0, to] — R, to > 0, be a continuous function. We say that u is an 
oscillatory function near the origin if there exists a sequence Sk such that 
Sk \ as — 7- oo, n(sfc) = and restrictions w|( Sfc+1 , Sfc ) intermittently change 
sign, k G N. 

Suppose f : / — ?• R, I — (0,1], is an oscillatory function near the origin, 
d G [1,2). We say that v is d-dimensional fractal oscillatory near the origin if 
dxm B G(v) = d and < M.* d {G{y)) < oo, where G{v) denotes the graph of 
function v. 

Let us define X : (0, l/t ] -)■ R by X(r) = x(l/r). We say that X(r) is 
oscillatory near the origin if x = x(t) is oscillatory near t = oo. We measure 
the rate of oscillatority of x(t) near t = oo by the rate of oscillatority of 
X(t) near r = 0. More precisely, the oscillatory dimension dim osc (x) (near 
t = oo) is defined as the box dimension of the graph of X{r) near r = 0. In 
Radunovic, Zubrinic and Zupanovic [TB] box dimension of unbounded sets has 
been studied. 

Assume now that x is of class C 1 . We say that x is a phase oscillatory function 
if the following stronger condition holds: the set T = {(x(t), x(t)) : t G [to, oo)} 
in the plane is a spiral converging to the origin. 

By a spiral here we mean the graph of a function r = f(ip), <p > (pi > 0, in 
polar coordinates, where 

f : [<pi, oo) (0, oo) is such that f(<p) — > as ip — > oo, 
< / is radially decreasing (ie, for any fixed tp > cpi (4) 
the function N 3 k i— > f(<p + 2A:tt) is decreasing), 

which is the definition from [T9] . 
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The phase dimension dim p h(x) of the function x(t) is defined as the box di- 
mension of the corresponding spiral T = {(x(t),x(t)) : t G [to? oo)}. 

We use a result for box dimension of graph G(X) of standard (a, (3)-chirps 
defined by 

X a Ar) = r a sm(r- 13 ). (5) 



For < a < f3 we have 



dim s G{X atP ) = 2 - (a + + 1), 



(6) 



and the same for X a ^{r) = t q cos(t~^), see Tricot [HJ p. 121]. Also we use 
a result for box dimension of spiral T defined by r = ip~ a , ip > y^o > 0, 
dim^r = 2/(1 + a) when < a < 1, see Tricot [HJ p. 121] and some 
generalizations from [19]. 

Oscillatory and phase dimensions are fractal dimensions, which are well known 
tool in the study of dynamics, see survey article [23] . 

For two real functions f(t) and g(t) of real variable we write f(t) ~ g{t), and 
say that functions are comparable as t — >■ (as t — > oo), if there exist positive 
constants C and D such that C f(t) < g(t) < D f(t) for all t sufficiently close 
to t = (for all t sufficiently large). For example, for a function F : U — > V 
with U,V C 1R 2 , V = F(U), the condition - F{t 2 )\ ~ |ti - t 2 | means 

that / is a bi-Lipschitz mapping, i.e., both F and F _1 are Lipschitzian. 

We say that function f is comparable of class k to power t~ a if / is class C k 
function, and f^\t) ~ t~ a ~ j as t — > oo, j = 0, 1, 2, . . . , k. 

Also, we write fit) ~ g{t) if f{t)/g[t) — )■ 1 as i — >■ oo, and say that function 
f is comparable of class k to power t~ a in the limit sense if / is class C k 
function, f(t) ~ t~ a and fV)(t) ~ (-l) 3 'a(a + l)(a + j - l)r Q -J' as t -> oo, 
j = l,2,...,fc. 

We write /(t) = 0(g(t)) as t — > (as t — > oo) if there exists positive constant 
C such that \f(t)\ < C\g{t)\. 

In the sequel we shall consider the functions of the form y = p(r) sin(g(r)) or 
V = P( r ) cos (q{ t ))- If P{ T ) — tOL i q{ t ) — T _/3 j a '{ T ) — t^ 13 ^ 1 as r — > then we 
say that y is an (a, (5)- chirp-like function. 
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3 Spiral trajectories in M 3 



In this section we describe solutions of equation ([T]) and trajectories of system 
with respect to box dimension, specifying a class of functions p. Let p(t) 
be comparable to power t~ a , a > 0, in the limit sense. Depending on a, we 
have rectifiable spirals with trivial box dimension equal to 1, or nonrectifiable 
spirals with nontrivial box dimension greater than 1. The box dimension will 
not exceed 2 even in IR 3 , because these spirals lie in a surface. Mapping spiral 
from the plane to the Lipschitzian surface does not affect the box dimension, 
see [2U1I2T] , while mapping to Holderian surface affects the box dimension. 

In order to explain fractal behavior of the system ^ we need a lemma dealing 
with a bi-Lipschitz map. The idea is to use result about box dimension of a 
class of planar spirals from [8, Theorem 4]. It is well known result from [2], 
that box dimension is preserved by bi-Lipschitz map. Putting together these 
two results we will obtain desired results about (|2]). For the sake of simplicity, 
we deal with trajectory T of the solution of the system defined by 



x(t) = p(t) smt 

y(t) =p'(t) sin t + p{t) cost (7) 

Trajectory Y has projection Y xy to (x, ?/)-plane which is a planar spiral satisfy- 
ing conditions of the mentioned |B1 Theorem 4]. In the following lemma we will 
prove that the mapping between planar spiral Y xy and spacial spiral Y is bi- 
Lipschitzian near the origin. We prove lemma using definition of bi-Lipschitz 
mapping. Interesting phenomenon appeared in spiral Y xy , defined in polar co- 
ordinates and generated by a chirp. The radius r{(p) is not decreasing function, 
there are some regions where r((p) increases causing some waves on the spiral. 
We introduced notion of wavy spiral in [8]. Also, the waves are found in the 
spiral generated by Bessel functions, and by generalized Bessel functions, de- 
pending on the parameters in the equation, see [9]. Furthermore, the surface 
containing the space spiral Y contains points with infinite derivative, showing 
some vertical regions. 

Lemma 1 There exists map B : M 2 x {0} — > M 3 such that B(x(t),y(t),0) = 
(x(t),y(t), z(t)) , where coordinate functions are given by and p(t) is com- 
parable of class 1 to t~ a , a > in the limit sense, as t — )■ oo ; and B is 
bi-Lipschitzian near the origin. 

Proof. It is clear that ^(r^) = Y, where Y xy is the projection of Y to (x,y)- 
plane, and Y is defined by parametrization (x(t),y(t), z(t)) from As from 
Theorem 4] it follows that Y xy does not have self intersections near the 
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origin, it is sufficient to prove that B\r is bi-Lipschitzian near the origin. 
That means we have to prove that there exist two positive constants Ki,K 2 
such that 



K 1 d((x(t 1 ),y(t 1 ),0),(x(t 2 ),y(t 2 ),0)) < (8) 
d{B{x{h), y(h),0), B(x(t 2 ),y(t 2 ), 0)) < 
K 2 d{{x{t 1 ),y{t 1 ),0),{x{t 2 ),y{t 2 ),0)), 

where d is Euclidian metrics and ti,t 2 > to, for t sufficiently large. Notice 
that, without loss of generality, t\ < t 2 . It is obvious that by Ki = 1 the left 
hand side inequality is satisfied. In order to prove right hand side inequality, 
first we prove 

(zih) - z{t 2 )f < C ((sfa) - x{t 2 )f + (yih) - y{t 2 )f) , (9) 



and then the right inequality will be satisfied. From planar case, [8, Theorem 
4], we know that 

v(t) = t + ^ + 0(r 1 ) = ~+t(i + o(r 1 ), t^oo (10) 

f(<f)-<f~ a , V^oo. (11) 



From [9, Proposition 1] and the fact that p(t) is comparable to t a in the 
sense of limit, it follows that there exists C\ G (0, 1), such that for every A<p, 
| < Ay? < 2tt+ |, holds 

for ip sufficiently large. Let 

Vl = <p(t 1 ) = n/2 + t 1 (l + 0(tf), (12) 
V>2 = <p(t 2 ) =7r/2 + t 2 (l + 0(t 2 - 2 ). 



Let us consider the case |y?2 — Vil < f • From P, Lemma 7], and (10), (12), we 
have 

y/ixfo) - x{t 2 )f + (yih) - y{t 2 )f > l -(y 2 - ^)fM > C 2 (t 2 - h)t 2 a . 
Hence, 



1 2^x2 _ ~ t\) 2 C 2 t 2 2 C 

l\ l 2 V\l 2 Lj 2 L 2 



(z(h) - z(t 2 )y = (-- -y = v \ 2 J J < (i3) 
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L/ 2 L l L 2 



< 



^ 2 ((^l)-^2)) 2 + (y(tl)-y(t 2 )) 2 ). 

For the case 2tt + | > \ip 2 — y?i| > f , we have 

V(^i)-x(t 2 )) 2 + (y(*i)-y(t2)) 2 > /M - fM = 

fM - f(<Pi + (^2 - > Ci(y? 2 - ^i)a^r Q_1 > 



Then again 



(14) 



(t 2 -ti) 2 

/2/ 2 

\2 i /„,/+ \ „,/+ \\2 



(z( tl ) - z(t 2 )) 2 = 



(xih) - x(t 2 )) 2 + (yih) - y(t 2 )Y 



< 



< 



r< 2 +4j-- 2 a- 2 

K 2 ((x( tl ) - x(t 2 )) 2 + ((y( tl ) - y(t 2 )) 2 ) . 
For the case \ip 2 — ipi\ > 2n + |, we define n : = 



V2-Vl-o 



2?T 



(15) 



Then we have 



'Ah) - x(t 2 )) 2 + (yih) - y(t 2 )) 2 > fM - fM 



n-1 



(f(<Pi + ^2tt) - /(^ + + 1)2tt)) + /(^ + n2vr) - /(p 2 ) > 



i=0 



n-l 



^ 27raC 1 (v?i + i2n)- a ~ 1 + -aCi(y?i + n2n)- a - 1 > 



i=0 



X>i + <27r)— 1 = ^ £ (2vr)— 1 (^ + <)— 1 > 
6 i=o 6 i=o Z7r 



r+n-l 



c$t l a 1 {t 2 — ti). 



fx 

2ir 



Furthermore 



(16) 



(z( tl ) - z{t 2 )f = 



(h-h) 2 

f 2 f 2 
L 1 L 2 



< 



(xih) - x(t 2 )) 2 + (yih) - y(t 2 )f 



2+ i+-2a-2 



C 5 t 1 t 1 



< 



(17) 
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K 2 {{x{h) - x{t 2 )f + {y{t x ) - y{t 2 )f) . 



From (13), Ubv, (17), the right hand side inequality d8h follows. □ 



Theorem 1 (Trajectory in IR 3 ) Let p(t) be a function comparable of class 3 
to power t~ a , a > 0, in the limit sense, as t — > oo. 

(i) Phase dimension of any solution of the equation ([!]) near the origin is equal 
to d\m ph (x) = J ^ a foraE (0, 1). 

(ii) Trajectory T of the system ^ near the origin has box dimension dims T = 
^ for a E (0,1). 

(hi) Trajectory T of the system for a > 1 is a rectifiable spiral and dims T = 
1. 



Proof. 

(i) For the sake of simplicity we take solution x(t) = p(t) sin t of the equation ([T]) 
and z — 1/t. Taking P(r) = p(t), t = 1/t, from [8, Proposition 1] we know 
that P(t) is comparable of class 3 to power r Q in the limit sense, as r — > 0. 
Spiral trajectory T of the system (|2| is defined by (J7|. Let T xy , defined by 
r = f((p), be the projection of T in (x, y)— plane. From [8j Theorem 4] we 
obtain that f(ip) ~ <p~ a , en G (0, 1) and dim B T xy = dim ph (x) = 

(ii) The map B : (x(t),y(t), 0) — > (x(t),y(t), z(t)) is a bi-Lipschitz map near the 
origin, see Lemma [I] It is clear that T = B(T xy ). The box dimension of a 
set is preserved under bi-Lipschitzian mappings, see Falconer [21 p. 44], so 
it follows dime T = j^—. 

(iii) The rectifiability of a trajectory of the system rt2J) for a > 1 can be proved 
directly. Indeed, from ^ and ([7]) for the length l(T) of the trajectory T we 
have 

oo 

/(r) = J ^/x 2 {t)+y2{t) + z 2 {t)dt < 
to 

oo 

J ^Jcft- 2a + Cit~ 2a + t-*dt 

to 

If -2a > -4 then 

oo 

J(r) <c J r a dt < oo, 

to 

and if —2a < —4 then 



J( r ) < C j r 2 dt < oo. 



to 



□ 
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Figure 1 System (k2J) for p(t) = t 4 , Lipschitz case 




case 



Figure 3 System Q for pit) = t 3 , Holder case 

Remark 1. The Bessel equation of order v has phase dimension equal to 4/3, 
for the proof see [9] . This is a consequence of a fact that the Bessel functions 
behave like chirp X(r) = r 1 / 2 sinl/r. 

Remark 2. The asymptotics of the first return map, also called the Poincare 
map near singularities of focus type, could be computed using the phase di- 
mension. The leading term for planar spiral from Theorem [l] (i) has power 
1 + 1, a > 0. 
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Theorem 2 (Projections) Suppose the function p(t) is comparable of class 3 
to power t~ a , a > 0, in the limit sense, as t — >■ oo. 

Then the projection G yz of a trajectory of the system ^ to (y, z)— plane is a 
(a, 1) -chirp-like function, and dime G yz = 3 - = ^ if a & (0, 1). 

Proof. Projection is 



, N . fi\ i /i\ i / ITTT ^7T7 A p( 1 ) \ 

y(z)=p - sin — hp- cos - = \ p (-) + p (-) sin — harctan — ttt)- 



For the functions P(z) = yV 2 (^) + P 2 {\) and Q(z) = ~ + arctan we 

have P(z) — z a , P'(z) — z a ~ l , Q(z) ~ z~ x , Q'(z) ~ z~ 2 as z — > 0. So y(z) is 
(a, l)-chirp-like function. To calculate the box dimension of G yz we apply [U 
Theorem 4]. □ 

Remark 3. In other words an oscillatory dimension of the solution of ([!]), 
under assumptions of previous theorem concerning pit), is equal to dim osc a; = 
3=2 if a £ (0,1). 

Natural question concerns rectifiability on the Holderian surface. For system 
([2]), a > 1 the Holderian surface did not affect the rectifiability. The following 
theorem gives a sufficient conditions for more general situation, a spiral lying 
into the Holderian surface z = g(r), g(r) ~ r 13 , (3 £ (0,1). Spiral is called 
Holder-focus spiral if it lies in the Holderian surface, and tend to the origin. 

Theorem 3 (Rectifiability in IR 3 ) Let f : [y?i,oo) — > (0, 00), ipi > 0, f((p) ~ 
f~ a > (/'(v 9 )! — CyT a ~ x , a > 1, r = f(<f) define a rectifiable spiral. Assume 
that g : (0, f((fi)) —> (0, 00) is a function of class C 1 such that 

g(r)~rP, \g'(r)\ < Dr^ 1 , /3 € (0, 1). (18) 



Let r be a Holder-focus spiral defined by r = f(tp), tp £ [<pi,oo), z = g(r), 
then T is rectifiable spiral. 

Proof. The corresponding parametrization of spiral T in Cartesian space co- 
ordinates is 



x = /(y?)cosy?, 

V = f(jp) Bin <p (19) 

z = g(f{<p))- 
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For the length l(T) of this spiral we have 

oo 

Z(T) = / sl&{<p)+i?{<p) + P{<p)dtp = 

oo 

/ ^/f 2 (<p) + r %) + g' 2 U^))P^)d v < 

oo 

C J \fp~ 2a + if- 2a ~ 2 + (f-^-^^-^dif = 

oo 

If -2a ft - 2 < -2a then 

J(r) < cy" yr a cfy> < oo, 

and if — 2af3 — 2 > —2a then 

oo 

1(F) <C J (p-^dy < oo. 
<pi 

□ 

Theorem 4 (New parametrization) Let x(t) = p(t) sin t, wherep(t) is a func- 
tion comparable of class 3 to power t~ a , a > in the limit sense, as t — > oo. 

By V we denote the spiral with parametrization (x(t),x(t),t~' y ), 7 > ; near 
the origin. Then 

(i) dim B r' = for a < 7 < I, 

(ii) r" is rectifiable and dims r" = 1 for a > 1, 
(hi) dim B r = 2 - for 7 < a < 1. 

Proof. Direct computation from the parametrization shows that for 7 > a, we 
have the Lipschitz case and for < 7 < a the case is Holderian. Taking into 
account that spiral comparable to r = ip~ a , a > 1 is rectifiable, using [2T| 
Theorem 7 and 9] or [201 Theorem 3.1. and Theorem 3.2.], and also Theorem 
[3] for rectifiable case, we obtain the claim. 

□ 

Remark 4. In the Theorem [4] the box dimension of the spiral V has been 
computed, and all values satisfy < dim^ r" < 2 — a for a G (0, 1). 
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4 Hopf bifurcation 



In this section we study the Hopf bifurcation of a class of systems with spiral 
solutions. Limit cycles are interesting object appearing in differential equa- 
tions. In particular, we consider a system having its linear part in Cartesian 
coordinates with a conjugate pair ±ui of pure imaginary eigenvalues with 
u > 0, and the third eigenvalue is equal to zero. The corresponding normal 
form in cylindrical coordinates is: 



f = a x rz + a 2 r 3 + a 3 rz 2 + 0(\r, z\) 4 

ip = ou + 0(\r,z\) 2 (20) 
i = b x r 2 + b 2 z 2 + b 3 r 2 z + hz 3 + 0(\r, z\) A , 

where a, and hi G R are coefficients of the system. Such systems and their 
bifurcations are treated in Guckenheimer- Holmes [31 Section 7.4]. The fold- 
Hopf bifurcation and cusp-Hopf bifurcation have been studied in Harlim and 
Langford [4J and the references therein. Here we make a note about box di- 



mension of a spiral trajectory of the simplified system (20) at the moment of 
Hopf bifurcation. 

In [151 . [22] we studied planar system consisting of first two equations from 



(21), and made fractal analysis of the Hopf bifurcation of the system. We 
proved that box dimension of a spiral trajectory becomes nontrivial at the 
moment of bifurcation. The Hopf bifurcation occurs with box dimension equal 
to 4/3, furthermore degenerate Hopf bifurcation or Hopf-Takens bifurcation 
occurs with the box dimension greater than 4/3. The more limit cycles have 
been related to larger box dimension. Analogous results have been showed for 
discrete systems in [5], and applied to continuous systems via Poincare map. 
On the other hand in [20] and [21] 3-dimensional spirals have been studied. 
Here we apply the results about the box dimension of 3-dimensional spirals to 



Hopf bifurcation appearing in (21) in cylindrical coordinates 



r = r(r 21 + E;=o a^ 2i )> 

= 1 (21) 

1 = b 2 z 2 + ... + b n z n . 



First two equations are standard normal form of codimension /, where the 
Hopf-Takens bifurcation occurs, see [17] . The third equation gives us the case 
where spiral trajectories lie on Lipschitzian or Holderian surface, depending on 
the first exponent. The Holderian surface has infinite derivative in the origin, 
geometrically it is a cusp. 
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We proved for the standard planar model that the Hopf bifurcation occurs 
with box dimension equal to 4/3 and the Hopf-Takens occurs with larger 
dimensions. Here we prove that on the Holderian surface 3-dimensional Hopf 
bifurcation occurs with the box dimension greater than 4/3. 



Theorem 5 (Hopf bifurcation) Let I = 1 in the system (21) and b p < be 
the first nonzero coefficient in the third equation and oq = then a trajectory 
T near the origin has: 



(i) if2<p<3 then 



(ii) if p > 4 then 



dim B r = -, 
3 



dim B r= 3 -— . (22) 
2 2p K J 

Proof. Using [191 Theorem 9] we get the solution of the first two equations of 



(21), r ~ cp 1 I 2 having dims = |, where T xy is orthogonal projection of 

— 2 

space trajectory T to (x,y) plane. From the third equation we get z ~ r?- 1 , 

so for 2 < p < 3 we get the Lipschitzian surface, while for p > 4 surface is 
Holderian. Applying [2H Theorem 7 (a)] we obtain dime T = 4/3 for 2 < p < 
3, because the box dimension is invariant for the Lipschitzian case. For the 
Holderian case we apply [211 Theorem 9 (a)], where a = 1/2 and = 2/ (p—1). 
So, we get dims r = § - ^. 

□ 



Remark 5. The box dimension of a trajectory at the moment of planar Hopf 
bifurcation is equal to 4/3, also for 3-dimensional case with spiral trajectory 
lying in the Lipschitzian surface. Situation is different for spiral trajectory 
contained in the Holderian surface, so the box dimension of a trajectory for 
the Hopf bifurcation tends to 3/2. Only one limit cycle could be produced, 
but dimension increases caused by the Holderian behavior near the origin. 



Notice that if we apply formula (22) obtained for the Holderian case, to the 
Lipschitzian case p = 3 we will get correct result 4/3. For / > 1 degenerate 
Hopf bifurcation or Hopf-Takens bifurcation appears, where / limit cycles could 
be born, and the box dimension is equal to dim^r = ^ 4 '~ 1 ^~ 2<+1 using the 
same arguments. 
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